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Abstract 

Causes of the unsatisfactory condition of the gravitational-wave experiments are discussed 
I and a new outlook at the detection of gravitational waves of astrophysical origin is proposed. 

It is shown that there are strong grounds for identifying the so-called giant pulses in the pul- 
' sar NP 0532 radiation with gravimagnetic shock waves (GMSW) excited in the neutron star 

, magnetosphere by sporadic gravitational radiation of this pulsar. 
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1 Introduction 



In the history of physics of the 20th century, I suppose, there is no such a grave experimental problem 
(except the controlled thermonuclear fusion problem) that, being solved for over thirty years by 
{2JQ' different research groups as the gravitational wave (GW) detection problem. Although much means 

are used to solve it, no sufficiently convincing possitive results have been obtained. What are the 
reasons for this situation? An error in the gravitational theory? The experimentalists' incapability? 
Are there realizable opportunities to detect gravitational radiation in the visible future? Is the GW 
detection problem worth studing? We will try to answer these questions in the present paper. As 
. any other radiation detection problem, this one also splits into two independent problems: (1) "GW 

sources" and (2) "GW detectors" . It is important not to forget to join both branches by solving a 
concrete experimental problem. 
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2 Sources of gravitational radiation 

2.1 Estimates of gravitational radiation power 

The average power of gravitational radiation from a source is calculated by the formula |9] 



Law — <tiktik>, (1) 
where G is the gravitational constant, c is the velocity of light, 

Uk=jp{...k-\8.kr-)dV (2) 

is the reduced quadrupole moment of the source; dots mean time derivatives. There are two param- 
eters of interest in the GW detection problems: the GW magnitude (i.e. the deviation from the flat 
metric hik = gik — Vik) and the GW frequency, uj. The GW energy flow is expressed in terms of these 
parameters by the formula 9 : 
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(Throughout the paper we use the metric signature (—1, —1, —1, +1)). From (1) - (3) follow estimat- 
ing formulae for the gravitational radiation power and magnitude: 



Law = ^f^; (4) 

Law = h , (5) 

where Eq is the energy of quadrupole ole oscillations of the source of the characteristic frequency a; 
{Lq = ijjEg is the quadrupole oscillation power), R is the distance from the source to an observer. 

2.2 Restrictions on GW magnitude 

In particular, a useful formula follows from ([4]) - ([5|): 

h E„ 
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where ho is the gravitational potential of the source of the total mass M: 



(6) 



ho - (7) 

According to ^ and (O), the ratio ofGW magnitude to the Newtonian gravitational potential of the 
source is of the order of the ratio of the quadrupole energy of the source oscillations to its complete 
energy at rest, Eq = Mc^. It is obvious that always Eg < Eq, and Eg/ Eq ^ 1 in a generic situations, 
therefore Eq. ([7]) gives an upper limit of the GW magnitude of a source, which is a good sobering 
factor by itself. Let us list some values for reference. For the Solar mass {Mq — 2 • 10^'^g) and a 
distance of 1 pc (3.26 light years — 3 • 10^^ cm)l, from Eq. ([7]) we obtaiii^^ 

/io-(pc/Af0)=4,8-lO-i^ 

For a mass of 1 kg at a distance of 1 m: 

ho ■ (m kg"^) = 7,4- 10^1^ 

Therefore for an A-bomb explosion with (AM/M ~ 10"'^), under the condition that the whole 
explosion energy turns into the quadrupole oscillations energy, at a distance of 1 m (!) from the 
epicentre, we get the GW magnitude h ~ lO"^""^. If one makes all atoms of a compact graser oscillate 
in the optic range (the radiation energy is about (hi/ ~ 1 ev, I ~ 1 m), we get the following maximum 
estimate for the GW magnitude on the graser end-wall: h ^ 10^^*. 

GW sources can be divided into two classes: (1) stable (quasistable) sources, which cannot be 
destroyed during the GW radiation process; (2) catastrophic sources, being destroyed in the GW 
radiation process. A graser represents a source of the first type, an A-bomb a source of the second 
type. For first type sources the quadrupole oscillation energy cannot exceed the binding energy of 
the source as a whole, unlike second-type ones, which are sources for one occasion. For example, 
close binaries and quadrupole oscillations of neutron stars are first-type astrophysical sources, and 
Supernovae are second-type ones. 



^The nearest stars' distance is about 1,3 pc. 



2 



As mentioned above, stable radiation sources are subject to the condition 

■^kin < ^b, (8) 

where E-^^^ is the inner kinetic energy of separate parts of the source, Ef, is their binding energy. 
Since it is always Eg ^ ^'kinj the condition ([8] takes the form 

Eq^Eb. (9) 

Hence the upper limit of GW magnitude from such sources can be obtained from the formula 
Eh 

h < ho^. (10) 

For astrophysical sources the binding energy is essentially that of gravitational attraction. Let AM be 
the part of the mass of an astrophysical object performing quadrupole oscillations. Its gravitational 
binding energy is 

^AM-M , , 

Eh < G J , (11) 

where / is the chracteristic size of the system. Thus for the upper limit of a GW magnitude from 
such a source Eqs. pH)) and pT|) give: 

where rg — 2GM/ is the gravitational radius of the radiating system. 
2.3 Radiation frequency 

Consider first a source of total mass M, which consists of two parts, so that the second part AM 
performs a free motion in the gravitational field of the system (rotation or free fall). Let w be a 
characteristic frequency of this procesfH. Equating the centrifugal and free fall accelerations, we 
obtain the well-known relation 

GM = uj^P, (13) 

which connects the characteristic size of the system with its characteristic frequency. 

Now let the gravitational attraction in the system be held by the forces of pressure (for stellar 
quadrupole oscillations). Equating these forces, we get the hydrostatistic balance condition 

, - , GM , , 

|VP| = p— , (14) 

where P is the pressure and p is the density. Using the known relationdP — v^^dp, where is the 
velocity of sound, from Eq. (IT4l) we obtain: 

Iv) w MG. (15) 



^Evidently the order of magnitude of this quantity coincides with the frequency of gravitational radiation from the 
system. 
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But Vf/l ~ Lo is the system proper oscillation frequency. Therefore for systems supported by the 
forces of pressure we return to the estimate of Eq. (13). 

Thus for stable astrophysical GW sources Eq. (13) has a universal nature if omega is understood 
as a characteristic frequency of the system oscillations. 

From the law (13) we can estimate the radiation characteristics of the collapsing objects, colliding 
stars and the like. It follows from this law that a maximum radiation frequency can be achieved for 
objects close to the gravitational collapse condition. In this case, by (12), the maximum magnitude 
of radiated GW is achieved (see [T]). For objects of masses of the order of the Solar mass (r^ —2, 96 
km) the maximum radiation frequency is 

(^max ^ c/rgf^ lO^sec"^ 



2.4 Close binary stellar systems 

For intense astrophysical sources of gravitational radiation, this radiation is the basic mechanism 
of quadrupole oscillation energy loss. Therefore, more rigorously, such sources should be called 
quasistable. The gravitational radiation power of a system of two orbiting gravitating masses ml 
and m2 is calculated from the known formula |10j 

_ dE _i2G^mlml{mi+m2) 

5^VS ' ^^^^ 

where r is the separation of the centres of mass. The energy balance leads to the mass approaching 
law [10] 

^ _ 64G^mim2(mi -I- TO2) 

Its integration yields a formula for the time t needed for the mass centres to approach to a distance 
of r from rp : 

^ = 192G3mim2(mi + ma) " ''o^' ^^^^ 

Further for simplicity we will study a pair of equal stars, setting mi — m2 = M , ro = 2Rq where 
Ro is the stellar radius, i.e. we will calculate the time until the catastrophic stellar collision, r (the 
lifetime). Then from ^ - (HH) we get: 

_ bc^jl^ - 16R^) 
~ 384G3M3 
and for I ^ i?o 

5 f I 



The gravitational radiation frequency of a binary system increases with time; the ratio of the fre- 
quency shift per period Aw to the radiation frequency u) is, by order of magnitude, 

--(f) ■ (20) 

Figs. 1 and 2 show the dependence of the gravitational radiation power of a binary system and 
the radiated GW magnitude on the distance between the stars. Fig. 3 shows the binary lifetime 
versus their separation for stars with the masses mi = TO2 = Mq. 
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Figure: 1. Gravitational radiation power of a binary Lg (erg/s) 
vs. distance between the stars I (km). 
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Figure: 2. GW magnitude h from a binary vs. its size I (km), at 
a distance of 1 kpc 

Let us estimate the probability oi GW detection from a close binary system in the Galaxy at 
given rotation period, assuming that the Galaxy age is of the order of 1 • 10^° years. Further, we 
assume that the average stellar number density in the Galaxy is of the order of 0, 120 stars/pc'^ [TS] . 
the Galactic volume is 300 kpc'^ [1^, then the number of stars in the Galaxy is about 0,35 • 10^^. 
Besides, we take into account that approximately half of the stars are in binary systems [17j. Then 
the possibility of existence of a binary with a prescribed lifetime r is proportional to the ratio r/t, 
where t is the age of the Galaxy. In the columns of Table 1 corresponding to system lifetimes 
smaller than 1 year, the detection probability of such systems in experiments lasting 1 year is shown. 
Evidently for such systems the probability of detection in a year-lasting experiment coincides with 
that for a binary having a lifetime of 1 year. The number of such systems in the Galaxy is estimated 
to be of the order of one. 
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Figure: 3. Lifetime of a binary r (seconds) vs. its size I (km) 

The presented data show that at the instant preceding the catastrophic colhsion, the gravitational 
radiation power from the binary is of the order of Supernova luminosity. Thus, a stellar collision in 
a close binary is an event whose scale is of the order of a Supernova explosion. As mentioned above, 
in the Galaxy the probability of GW detection from a binary with a lifetime of the order ofl year is 
close to one. This means that catastrophic phenomena with energy release of the order of 1 • 10^'* 
erg/s should happen once a year. However, in reality such phenomena happen once in 40 to 80 years 
in the Galaxy [IQj. A reason for such a discrepancy is that in a close binary with a lifetime of the 
order ofl year the stars' separation is about 12000 km. 

Therefore for such a system to exist it is necessary that both stars be at least white dwarfs. But 
in this case even in a stellar collision the energy released is 4 orders of magnituded smaller than 
that of a Supernova explosion. Thus for a catastrophic collision of this scale it is at least necessary 
that one of the components be a neutron star, while the second one is a white dwarf. The existence 
probability of such systems in the Galaxy is much smaller. Note that it is difficult to understand 
the experimental programmes intended for registration of GW from the binaries with periods of the 
order of a few seconds. According to Table 1, their lifetime does not exceed 5 years, and in this case 
it would be more reasonalble to wait these 5 years and to detect the gravitational radiation from a 
catastrophic collision: its power is higher by at least 13 orders and the GW magnitude is greater 
by 3 orders (!), as follows from Table 1. However, at least in the last 10 years nobody detected 
catastrophic events on such a scale at distances smaller than 15 kpc. 

Since the pulsars are identified with Supernovae remnants, the average frequency of Supernova 
bursts may be estimated from the data on pulsars spreading in the Solar neighbourhood. Thus, 
at distances within about 1 kpc, on the whole, about 20 pulsars are observed. Table 2 shows the 
data on the pulsars nearest to the Solar system3. As follows from this Table, almost all the pulsars 
are younger than 108 years. Therefore, it may be stated that the observed pulsars are remnants of 
Supernovae which exploded in the last hundred million years. It gives 1 flash per 50 years, coinciding 
with the estimate of Ref . [5] . It seems likely to be also close to the average frequency of catastrophic 
collisions in close binaries. 
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Table 1: Characteristics of gravitational radiation from a close binary* 



I 


5(4) 


2(4) 


1(4) 


7(3) 


4(3) 


1(3) 


320 


100 


40 


20 


T 


8,2(10) 


2,1(9) 


1,3(8) 


3,15(7) 


3,3(6) 


1,3(4) 


129 


1,30 


0,031 





LU 


0,046 


0,18 


0,52 


0,88 


2,05 


16,4 


92,6 


518 


2047 


5790 


T 


136 


34,3 


12,1 


7,10 


3,07 


0,38 


0,068 


0,012 


0,003 


0,001 




1,1(38) 


1,1(40) 


3,4(41) 


2,0(42) 


3,3(43) 


3,4(46) 


1,1(49) 


3,4(51) 


3,3(53) 


1,1(55) 


h 


7,5(-21) 


l,9(-20) 


3,7(-20) 


5,4(-20) 


9,4(-20) 


3,7(-19) 


1,2(-18) 


3,7(-18) 


9,4(-18) 


1,9(-17) 


V 


9,7(-14) 


9,5(-12) 


3,0(-10) 


l,8(-9) 


3,0(-8) 


3,0(-5) 


9,8(-3) 


3,0 


296 


9463 


N 


2343 


60 


3,7 


1 


1 


1 


1 


1 


1 


1 




l,5(-20) 


4,4(-21) 


1,7(-21) 


3,7(-21) 


6,5(-21) 


2,6(-20) 


8,3(-20) 


2,6(-19) 


6,5(-19) 


1,3(-18) 




3,8(-13) 


3,2(-12) 


1,6(-11) 


9,0(-12) 


1,5(-10) 


l,5(-7) 


4,7(-5) 


l,5(-2) 


1,4 


46 



*^ Here and henceforth, the figures in parantheses indicate the order of magnitude (a(b) = a- 10* J; the 

quantities h and V (the gravitational radiation Bow density) are calculated at a distance of 1 kpc from the 
binary; T = 2-k/uj is the GW period; N is the expected number of binaries in the Galaxy; R =< R > is the 
expected distance to the binary in kpc; /i® is the expected GW magnitude on the Earth; P. is the expected 
gravitational radiation How density on the Earth. The quantities T and t are given in seconds, ui in sec~^, I 
in km, Lg in erg/sec, V and in W/cri? . 



Table 2: Data on pulsars located at distances smaller or of the order of 1 kpc from the Sun 



No 


Pulsar 
(name) 


Distance 
(kpc) 


P 

(sec) 


P/P 
(years 


1 


MP 0031 


0,21 


0,94 


7,1(7) 


2 


MP 0450 


0,33 


0,55 




3 


NP 0532 


2,0 


0,033 


2,5(3) 


4 


MP 0628-28 


0,170 


1,24 


1,6(7) 


5 


CP 0809 


0,19 


1,29 


2,5(8) 


6 


AP 0823+26 0,38 


0,53 


1,0(7) 




7 


PSR 0833-45 


0,5 


0,089 


2,3(4) 


8 


CP 0834 


0,43 


1,27 


5,9(6) 


9 


PP 0943 


0,30 


1,098 




10 


CP 0950 


0,10 


0,253 


3,5(7) 


11 


CP 1133 


0,16 


1,188 


1,0(7) 


12 


AP 1237-)-25 


0,20 


1,382 


4,6(7) 


13 


PSR 1451-68 


0,40 


0,263 i2,8(6) 




14 


HP 1508 


0,26 


0,740 


4,7(6) 


15 


CP 1919 


0,42 


1,337 


3,2(7) 


16 


PSR 1929-1-10 


0,27 


0,227 


6,2(6) 


17 


JP 1933-1-16 


3 


0,359 


1,9(6) 


18 


AP 2016+28 


0,47 


0,558 


1,2(8) 


19 


PSR 2045 


0,38 


1,962 


5,7(6) 
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Thus we have to deal at best with system Unear sizes of the order of 110000 20000km. It gives: 
T ~ 10-h 40sec, /le ~ (1 5) • IQ-^i, Pgj ~ 10"" ^ 3 • IQ-i^W/cm^. GW with such parameters can 
hardly be detected in the coming decades. In this situation it only remains to hope for a case, rare 
and simultineously dangerous for the Earth, of a Supernova burst or a catastrophic end of a close 
binary. 

2.5 Neutron star oscillations 

There is, however, one more class of stable astrophysical GW sources — quadrupole oscillations of 
neutron stars. Table 3 shows the calculated parameters of gravitational radiation from neutron stars 
1^ (columns 1 7). According to this table, the following characteristics of radiation are to be 
expected from these sources: ~ (0, 3 1) • 10'*s~^; ^ 10^^^, V@ ^ 3 • lO^^^W/cm^ by the energy 
ofstellar quadrupole oscillations of the order of 10'^*-^ 10^^ erg and the whole gravitational luminosity 
of Lg ~ 4 • 10639 erg/sec. Note that due to the smallness of the GW magnitude from this source 
and the absence of a mechanism able to support the excitation of a necessary quadrupole moment 
during a sufficiently long time, quadrupole oscillations of neutron stars have not been considered as 
a competitive GW source. 

3 GMSW and GW detection 

The cause of the unsatisfactory condition in the GW detection problem is, in the author's opinion, 
the originally chosen erroneous way ofits solution — the programme ofcreating GW detectors. Direct 
GW detection can be realized either due to their tidal effect on a nonrelativistic (solid-state) detector, 
or due to their relativistic effect on a detector having a relativistic component (a laser ray). In both 
cases the GW effect on a detector (test-body displacement or laser ray deviation) is propotional to 
the GW magnitude. And the expected GW magnitudes from astrophysical sources are extremely 
small (see, for example, [9]). 

The existing GW detection programmes are generally meant for astrophysics sources of two types: 
1. Supernovae; 2 close binaries. In the first case one may expect GW magnitudes about 10~^^-;-10~^* 
with the radiation in a wide frequency range with the characteristic frequency of the order of 10^ 
sec""'^, in the second case magnitudes about 10^^° 10^^^ with a fixed frequency in the range of 
0, 1 10 sec~^. Due to the very small expected GW magnitudes on the Earth, the experimental 
programmes intended for direct GW detection inevitably come across the problem of noise of external 
thermal and quantum nature. This struggle is already in its third decade and requires the creation 
of high-precision deeply cooled detectors. 

On the other hand, it is well-known that even such weak-magnitude GW carry rather a high 
energy: in the above examples, this energy is of the order ofl W/cm2 in the first case and about 
10-13 _^ 10~" W/cm in the second case. Electromagnetic signal detection on such a power level has 
no problems. Therefore, the GW detection problem should be solved in a different way: by looking 
for specific electromagnetic signals from GW effect upon matter in those regions of the Galaxy where 
the gravitational radiation intensity is high. Setting the problem in such a way, we should above all 
study the GW effect on plasma-like media. The corresponding studies, carried in the eighties mainly 
in the Kazan school of gravitation, revealed a number of specific electromagnetic reactions of plasma 
to GW. In Refs. [2] — [5] the effect of plane GW (PGW) on plasma-like media was investigated by 
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the methods of relativistic kinetic theory in the approximation of neghgible back reaction of matter 
on the PGW: 



(87rG/c2)e<a;2. (21) 

where uj is the GW characteristic frequency, e is the matter energy density. These papers have 
revealed a number of phenomena of interest, consisting in induction of longitudinal electric oscilla- 
tions in the plasma by PGW. In spite of the strictness of the results obtained in [5]-[S], the effects 
discovered have very little to do with the real problem of GW detection. Moreover, the above cal- 
culations show a lack of any prospects for GW detectors based on dynamical excitation of electric 
oscillations by gravitational radiation. There are two reasons for that: the smallness of the ratio 
m^G/e^ = 10^^'^ and the small relativistic factor {v'^)/c^ of standard plasmalike systems. The GW 
energy conversion coefficient to plasma oscillations is directly proportional to a product of these 
factors. 

However, the situation may change radically if strong electric or magnetic fields are present in 
the plasma. In Ref. [B], where the induction of surface currents on a metal- vacuum interface by a 
PGW was studied, it was shown that the values of currents thus induced can be of experimental 
interest. In fTl, on the basis of relativistic kinetic equations, a set of magneto- hydrodynamics (MHD) 
equations was obtained, which described the motion of coUisionless magnetoactive plasma against the 
background of a PGW of arbitrary magnitude in drift approximation and it was shown that, provided 
the propagation of the PGW is transversal, there arises a plasma drift in the PGW propagation 
direction. 

In Ref. [1] an exact solution of the relativistic MHD equations in the PGW background of 
arbitrary magnitude was obtained and, on its basis, a fundamentally new class of sufficiently nonlinear 
threshold effects was discovered, named GMSW ( "jimmysway" ) - gravimagnetic shock waves. 

3.1 GMSW 

The PGW metric of the polarisation e-|- is described by the expression [5]: 

= 2dudv - L^[e^^{dx^)^ + e-^^{dx')% (22) 

where f3{u) is an arbitrary function (the PGW magnitude), while L{u) (the PGW background factor) 
obeys the ordinary second-order differential equation 

L" + LI3^ = 0; (23) 

u = l/\/2{t — x^) is the retarted time and v = l/^/2{t + x^) is the advanced time. Let there be no 
PGW at (li «C 0): 

while the plasma be homogeneous and at rest: 

^''HL^o = «n")l„^o = l/%/2; vX = v\ = 0; 

(25) 

£(")l«<o; p(w)L<o = po 
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{p = P{^) is the plasma pressure, is its dynamic velocity vector) and a homogeneous magnetic 
field is directed in the plane: 

-^i('")lu<o = Hncosi}; iJ2(w)|„<o = HnsmV,; 

(26) 

where ft is the angle between the axis Ox^ (the PGW propagation direction) and the magnetic field 
H direction. The conditions (26) correspond to the vector potential : 

Ay = A^ = 0; 

As ^ Hq{x^ Sinn ~ x'^ cosn); (m ^ 0). (27) 

The exact solution of the relativistic MHD equations against the metrics background ((22]) . The 
exact solution of the relativistic MHD equations against the metrics background (1^^ obtained in [T] 
satisfies the initial conditions (24) - (26) and is determined by the governing function :) obtained in 
[I] satisfies the initial conditions (PS)) - (P7)) and is determined by the governing function: 

A{u)''J l-a\e^^ -1), (28) 
where a is a dimensionless parameter: 
2 sin^ n 

" ^ -r-7 29 

47r(£o +Po) 

This solution contains a physical singularity on the hypersurface S : u = m*: 

A(u,) = 1 - a2(M,)(e2^("*) - 1) = 0, (30) 

where the plasma and magnetic field energy densities tend to infinity and the dynamic velocity of the 
plasma as a whole tends to the velocity of light in the PGW propagation direction. In this case the 
ratio of the magnetic field energy density and the plasma energy tends to infinity. This singularity is a 
gravimagnetic shock wave (GMSW, [I4), spreading in the PGW propagation direction at a subluminal 
velocity. According to Eq. ([30)) . necessary conditions for the occurrence of the singularity are 



/3(u) > 0; (31) 

> 1. (32) 

An extrimely important fact is that a singular state is even possible in a weak PGW (|/3| <C 1) 

under the condition that the plasma is highly magnetized (a^ 3> 1); in this case the singularity 
condition arises according to (30) on the hypersurfaces u = u*: 

/3(u,) = l/(2a2). (33) 
In particular, for a barotropic equation of state (p — fee, ^ k < 1) 

e = £oA-'+"; (34) 

= -^L-'Ai+t; (35) 
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^ = A-2 rA-- + (A - l)2L-2e-2/3 cot^ fil ; (36) 



H^ = ^ = (cos^ ^1 + L^A-^e^'^ sin^ O) , (37) 



A2 

where 



and 



9k 

A = L\u)A{u), u=^>0, 
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is the electromagnetic field invariant, (squared magnetic field strength in the frame of reference 
comoving with the plasma) . 

It follows from (34) - (37) that if ^ > 0, the plasma moves in the GW propagation direction 

{v^ = l/\/2(fM — Wr) > 0) and if /3 < 0, in the opposite direction. The effect is maximum in the 
PGW propagation direction, that is, perpendicular to that of the original magnetic field strength, 
and vanishes in the direction parallel to the magnetic field strength. 

In the case of strictly transversal PGW propagation ($1 = 7r/2), in the direction Ox^ the plasma 
drift vanishes, and the component of the plasma physical 3- velocity in the Ox^ direction is 

1 Vu-Vy 1 - A^A" 

V =c = c- . „ , . (38) 

The components of the total (including the magnetic field) EMT of the magnetoactive plasma, 
(a = 1,4) have a hydromagnetic structure: 

T^ = {£ + P)v'^Vk-P5t, (39) 

where 

eH = PH = ^; £ = e + eH\ P = p + Ph, (40) 

OTT 

where P and E are the total pressure and energy density of the magnetoactive plasma. There arises 
an energy flow in the plasma in the direction Ox^ : 

= ^-^±^{A-^A^- - 1){AA^- + a^e^^). (41) 

The parameter v takes in these formulae the following values in the two extreme cases: 

0; fc = 0; 
1; A; = 1/3. 



(42) 



For a weak GW 

|/3(u)|«l; L\u) = l + 0{p^)^l, (43) 

the expressions (37), (41) and (38) take the form 

_ 1 - A™ 
c ~ 1 + A™ ' 
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= J(eo+Po)(l + a')(A-"-l); 



(44) 



1 



(45) 





where the coefRcients m and n take integer values for nonrelativistic {k = 0) and uhrarelativistic 
(fc = 1/3) plasma: 



4 GW energy transmission to plasma: a half-self-consistent 
solution 

4.1 Total momentum conservation 

Since on the singular hypersurface (30) A(m) = 0, the energy densities of the plasma and the magnetic 
field tend to infinity, and the velocity of the plasma as a whole tends to the speed of light, the total 
energy of the magnetohydrodynamic shock wave and its flow in the GW propagation direction tend 
to infinity. The singular state emerging in the plasma due to the PGW violates the basic assumption 
(21) of the weaknesss of GW interaction with the plasma. In a more complete self-consistent problem 
including gravitation, the back reaction of the shock wave upon the PGW should lead to PGW energy 
loss and its magnitude damping up to the values 



Thus a GMSW is an effective mechanism of a gravitational wave energy pumping over into plasma 
[T]. A rigorous solution of the problem of PGW energy transformation into the shock wave energy is 
only possible by studying the self-consistent set of the Einstein equations and the MHD equations. 

Ref. [T] suggested a semiquantitative solution of this problem on the basis ofa simple model 
of energo-ballance. Due to its extreme importance, we do not restrict ourselves to [T| and return 
to a more complete study of the problem of energy transmission from a GW to magnetoactive 
plasma. However, instead of solving the Einstein equations, we make use of their consequence, the 
conservation law of the total momentum of the system "plasma -I- gravitational waves" . Clearly, 
this model is only approximate and cannot replace a rigorous solution to the Einstein equations. 
According to IQi , an arbitrary gravitational field provides the conservation of the system's total 
momentum 



where f"^ is the energy-momentum pseudotensor of the gravitational field and the integration covers 
the whole 3-dimensional space. Let us take into account that the above solution is plane-symmetric 
and only depends on the retarded "time" u. Consequently the integration over the "plane" (x2,x3) 
in (|48|) reduces to simply multiplying by an infinite 2-dimensional area. Dividing both sides of (|48|) 



k = Q; m 
k — 1/3; m 



= 2; n = A- 
= 3; n = 2. 



(46) 



max|^| < l/2a^. 



(47) 




(48) 
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by this area and bearing in mind that with D, — tt/2 among the 3-dimensional flows only is 
nonzero, we obtain the conservation law of the surface density of the momentum P^: 



Ph^l I {-9){T'^ + t'^)x = Const. (49) 

Let the right semispace a; > be filled with magnetoactive plasma and the left one x < with 
matter which does not interact with a weak GW. Let further the whole gravitational momentum be 
concentrated in the interval u £ [0,u/], where tf — \/2uf is the gravitational pulse duration. Since 
the integral in Eq. (49) is conserved all the time, let us consider it at to < 0, when the GW has 
not yet reached the magnetoactive plasma, and tf > t > 0, when the GW has reached the plasma. 
Taking into account that the vacuum solution depends only on the retarded time, we get for the 
integral in Eq. (49): 

Uf Uf 

tl'= I {T^^ + t^^)du+ I tl\ (50) 

t/V2 

where ij^ = t^^{(3oiu))', — i^''(/3(u)), Poiu) is the vacuum magnitude of the PGW, (3{u) is the 
PGW magnitude with allowance for interaction with the plasma. Transferring one of the integrals 
to the left-hand side of Eq. (50), we arrive at the relation 

u 

j{T^^ + t^^)du, (51) 


where the variable u = t/^/2 > can now take any positive values. 

A similar law may be written for the plasma total energy; in this case instead of Eq. (51) we 
obtain: 



u u 





where Sq is the total energy density of the unperturbed plasma. 
4.2 Local analysis of the conservation law 

Since the relation (51) must be valid at any values of the variable u, the corresponding local relation 
should hold: 

T^'{/3)+t^\P)^t^\(3o), (52) 

i.e. a local conservation law of the energy flow density should hold, as was assumed in Ref. 1 . It 
should be pointed out that the local conservation law (52) is a direct consequence of the solution 
stationarity, i.e. the solution dependence on the retarted time u = {ct — x)/y/2- There are two factors 
preventing the solution in a rigorous model from being stationary: (1) PGW interaction with the 
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plasma; (2) the boundary conditions on the surface x = 0. In accordance with the approximation 
(21), we introduce a small dimensionless parameter x [I]: 

2 _ TTGjeo + po){l + a^) 
where uj is the characteristic GW frequency. 

The approximation (21) is equivalent to the condition 

X' « 1- (54) 

Under the condition (54) the GW velocity tends to that of light, thus providing the required solution 
stationarity even in inhomogeneous plasma [24 . Let /3o = Const > be a maximum value of the 
PGW vacuum magnitude, Let us introduce one more dimensionless parameter, the first GMSW 
parameter 

= xV/3o - ^I^GW (55) 

where 

Thus, the parameter ^ is of the order oft he ratio of the total magnetoactive plasma energy to the 
vacuum GW energy. 

Making use of the solution of the MHD equations in the case of strict transversal PGW propa- 
gation (il = P«/2) as well as the expression of the total plasma EMT (44) and that for the energy 
flow of a weak PGW (3), we reduce Eq. (52) to the form 

/32+F(/3)=/32, (56) 
where 

F(/3)=x'[A-"(/?)-l] (^. = 1) (57) 

is a function of /3; /3 is now a derivative in the dimensionless "time" s = ^f^bjujc. 
Introduce the relative PGW magnitude: 

q = /3//3o; = /?*//3o. (58) 
Then Eq. (56) may be rewritten in the form 

q^ + V[q)^ql. (59) 
where 

F(g)=e'[(l-Tg)-"-l] (60) 
and a new dimensionless parameter has been introduced: 

T = 2c?fil (61) 
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(the second GMSW parameter). The total energy conservation law leads roughly to the same result. 
Eq. (59) may be treated as an equation with respect to the variable q . On the other hand, (59) 
completely coincides in its form with the energy conservation law of a 1-dimensional mechanical 
system described by the canonical variables {g(s),g(s)} [23], where V{q) is the potential, is its 
kinetic energy and q^ = Eq is its total energy. Fig. 4. shows the qualitative form of the potential 
V{q). 



Figure: 4. Potential V{q) ofEq.(59) 

Two points on the potential curve (A and B) correspond to any positive value of Eq. These points 
are the system trajectory turning points. No real system states exist under the potential curve V{q). 
At the point 

g = 9c - T-i; (/? = l/2a') 



Viqc) 



oo. 



To analyze the system behaviour, let us suppose that the moment s 
edge of the GW, while 



/3* « (3q sin s : 



sin(s). 



(62) 
(63) 

corresponds to the front 
(64) 



Thus, provided the initial conditions (24) are satisfied, the system always starts from the point 5*0 
along the line (AB) towards A (for (3 > 0). Since this is a turning point, the maximum accessible 
value of the variable q in the system is q{A). This is the smallest root q- = (/(x, T,i?o) of the 
algebraic equation 



V{q) -Eo=0. 

The maximum attainable PGW magnitude in the system, /3„ 

Pmax — q-f^Q- 



(65) 



(66) 



Thus q_ coincides in its sense with the "PGW magnitude damping factor" 7 introduced Ref. [T]. 
Solving Eq. (j65p , we obtain the required root q^ : 



9- = Y 



1 + 1 



(67) 
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From it follows that always 

q- sS (68) 
and also, as i?o 

- ^ 0. (69) 

With increasing Eq this magnitude grows and for i?o — > oo it reaches the value {q = qc): 

After the turning point the GW magnitude diminishies, reaching negative values. For s — > +oo 

s -> +00 f3' f3'^ = Const < 0; /3 ^ -oo; 

the metric (22) degenerates (322 — >■ 0,(733 — —00); the only nonzero components of the curvature 
tensor take the following form due to the Einstein equations (see (23)): 

i?„2u2 = {L^y P'oo eM^P'oou) ^ -0; 

Ru3u3 = {L^y P'oo exp(2/?;,w) ^ +00. 

Thus, as s — ^ +00, a true singularity is formed in the system. It is easily verified that in this case 
— )■ 0, e — >■ 0, — — c. The plasma in the final state moves to meet the original GW direction, 
this reverse of the plasma needs a more detailed self-consistent analysis. 

4.3 Numerical analysis ofGMSW 

Let us pass to a more detailed study on the selfconsistent motion of the system. From Eq. (59) we 
obtain the differential equation 



dq/ds^±yM^V{q) (71) 

where the plus sign is chosen before and the minus sign after the turning point g_. It is helpful to 
solve and analyze Eq. (71) using the new dimensionless variables: A(/3) and 

S = Ts = V2Tuju (72) 

Substituting into (71), for example, q*{s) — sins, we reduce it to the form 



dA/dS = Tv/cos2 S/T~V{A), (73) 

moreover, the initial condition is to be fulfilled: 

A(0) = 1. (74) 

Figs. 5-11 show some results of numerical integration of Eq. (73) with the initial condition 
(74). An analysis of the formulae describing GMSW and numerical calculations make it possible to 
discover a number of general laws of the GMSW excitation process in homogeneous and isotropic 
plasma under the condition that the PGW propagation is strictly transversal: 
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1. A GMSW is completely described by three nonnegative dimensionless parameters: the param- 
eter k in the plasma equation of state, the first (f^) and the second (T) GMSW parameters. 

2. Necessary conditions for GMSW excitation are (31) and (32): 

T ^ 1. (75) 

3. The only criterion of strong GW absorption is, according to (69), a large value of the second 
GMSW parameter (T): 

T>1. (76) 

4. Under these conditions a maximum response of the plasma to GV^ is achieved when the values 
of the first GMSW parameter are small: 

< 1. (77) 

5. The plasma response to GW is a single pulse, and the shock wave stage is always replaced by a 
reverse stage, when the plasma turns back. Simultaneously its density, pressure and magnetic 
field strength fall off. 

6. The ultrarelativistic (fc — 1/3) plasma response is much greater (by approximately 2 orders) 
than that of a plasma with the nonrelativistic equation of state (fc = 0), and in ultrarelativistic 
plasma the pulse duration is also slightly greater. 

7. The profiles of the plasma response at sufficiently large values of the second GMSW parameter 
(T ^ 5) actually coincide on the S scale. This means that in the conventional time scale t the 
pulse duration is inversly proportional to the second parameter, or more precisely 

^^-2^ = 4^1' 
where T is the GW period. 

8. With T < 5 the response magnitude rapidly decreases and at T ^ 1 becomes smaller by an 
order of magnitude. In this case a maximum pulse duration is achieved: 

Ts^T/4 = 7r/(2a;). (79) 

9. A decrease in the first GMSW parameter causes a rapid increases in the response (roughly 
propotional to l/'C^); simultaneously increases the pulse duration approximately by a factor of 
2. 

10. A maximum response is achieved at the instant ~ 1. 



^Here and further on, speaking of a maximum response of the plasma, we mean its energy characteristics: the 
plasma energy flow density and the magnetic field energy density. 
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11. Under the optimal GMSW conditions (76) and (77), the total surface density of the magnetic 
field energy, transported in the pulse, fs, is f or ultrarelativistic plasma of the order of 

£^ ~ cifoVK'Tf)- (80) 

12. The shock wave energy is taken from the GW energy, so under the conditions (76) and (77) 
the GMSW is an effective mechanism of gravitational waves energy transformation into other 
forms of energy. 




Figure: 6. Physical velocity of plasma in GW Geld,v^{S)/c: 1 — - 
nonrelativistic, 2 — - ultrarelativistic equation of state; T = 10; 
^2 =0,01. 
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Figure: 7. Physical velocity of plasma in GW field, v^{S)/c: 1 - 
r = 10; 2 - T = 100. Everywhere = 0, 1. The lines practically 
coincide. 
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A I 1 > -I 

Figure: 9. Magnetic field strength Ig H'^{S)/H^ in GW field, 1 
-^k = 0,2^k = 1/3; T = 10; = o,01.. 




Figure: 10. Magnetic field strength in GW field for ultrarelativis- 
tic plasma, H'^{S)/H^: I ^ T = 10; 2 ^ T = 100. Everywhere 
= 0, 1. The lines practically coincide. 
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11 11", I . J 

Figure: 11. Magnetic Geld strength in GW field for ultrarelativis- 
tic plasma, lgH^{S)/H^: 1 — = 1; 2 — = 0,1; 3 ^ 
^2 = 0, 01. Everywhere T = 10. 
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Figure: 12. Magnetic field strength in GW held for ultrarelativis- 

tic plasma, H'^{S)/Hl 1 ^ T = 1; 2 T = 10; 3 

T = 100. Everywhere = 0,01. The lines 2 and 3 practically 
coincide. 



5 GMSW in neutron star magnetospheres 

5.1 GSMW parameters in neutron star magnetospheres 

In [1] it was shown that in the magnetospheres ofneutron stars performing quadrupole oscillations, 
large values of the second GMSW parameter are realized. That is, the necessary condition for GMSW 
excitation (76) is fulfilled. Let us study this problem in more detail. The electron number density in 
a pulsar magnetosphere, ne{r), which is necessary for calculating the parameter T, may be obtained 
by dimensional estimation from the Maxwell equations [12]: 

ne(r) ~ iJ(r)/(47rer). (81) 
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Further, as is known from (see [H!), the pulsar period slowing-down rate to is connected with the 
pulsar parameters as follows: 

where R is the neutron star radius, M is its mass, H is the magnetic field strength at the stellar 
surface, P is the rotation period. This formula gives for the pulsar NP 0532 NP 0532 77 w 5 • 10^^ 
G (see Table 3) according to the known slowing-down rate of this pulsar. Actually, as pointed out in 
Ref. [TT] , the magnetic field strength at the surface of NP 0532 is somewhat smaller than the value 
obtained on the basis of (82), and is of the order of 10^^ G. Further on we use this value. In the 
figures shown below (unless specially indicated) the following values of the parameters are adopted: 
R = 1, 2- 10^ cm, Po{R) — 10~^ and the magnetic field in the magnetosphere is assumed to be dipole: 
H{r) - (i?/r)3. 



Table 3. GMSW in a neutron star magnetosphere 



rJR 


SM 






Em/AR 




R 


H{R) 






T 


Em 


0,057 


0,405 


1,197 


13,0 


7,8(50) 


1,2(50) 
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7,7(11) 


l,0(-4) 


5,2(11) 


1,1(5) 


8,4(42) 


0,159 


0,677 


0,699 


1,7 


5,7(52) 


7,0(52) 
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2,7(12) 


4,2(-6) 


4,7(11) 


2,4(4) 


1,1(42) 


0,240 


0,682 


0,311 


0,2 


2,8(52) 


2,9(53) 
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6,1(12) 


3,3(-6) 


7,5(11) 


4,4(4) 


3,0(41) 


0,580 


1,954 


0,378 


0,2 


1,7(54) 


1,6(55) 


10 


7,4(12) 


4,2(-7) 


8,0(11) 


4,2(4) 


3,1(41) 


0,434 


1,670 


0,349 


0,2 


5,0(53) 


5,0(54) 
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5,2(12) 


9,5(-7) 


6,3(11) 


4,1(4) 


4,7(41) 



* Comments to Table 3. The data placed in columns 1 - 7, 9 and 12 are taken from the book [9]. 
5M = M /Mq is the neutron star mass related to the Solar mass; R is the star radius in km; T„ is 
the neutron star eigen-oscillation period in the basic quadrupole mode (in milliseconds); r„ is the 
oscillation damping time (in seconds); Ai? = {{SR/R)"^) is the root-mean-square relative magnitude 
of the neutron star oscillations; Em is the oscillations kinetic energy in erg; Lg is the the star's 
gravitational luminosity in erg/sec. The GW magnitude value at the neutron star surface, l3o{R), is 
assumed to be equal to 10~®. H{R) is the magnetic field strength in gauss (G). The data placed in 
columns 8-^11 are calculated using Eqs. (3), (82) and (81) for the observed Crab pulsar (NP 0532) 
parameters; P = 0,033 s, to — 2500 years. The data placed in the last line of the table (columns 2 
7 and 9) are obtained by extrapolation of the values from the book |9j . The values of the parameters 
and T are given for the magnetosphere near the stellar surface. 



If the magnetic field of a neutron star is described as that of a dipole, then the geographic angle 
6 (counted from the magnetic equator) will be connected with the above angle n by the relation 
= 7r/2 — 0. Therefore the GMSW excitation condition depends on the angle &: 



sin^ e < 1 - ^ l,^, - 1 - T"^ 
2al\f3\ 



Thus, in the magnetosphere of a neutron star (or a Supernova) a GMSW can be excited in the 
vicinity of the magnetic equator, similarly to pulsars, with a knife radiation pattern. In this region, 
as was demonstrated by the above examples, the gravitational radiation can be absorbed almost 
completely by shock wave excitation. Fig. 13 shows the radial dependence of the GMSW parameters 
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magnetic equator plane (O = 0) of a neutron star magnetosphere with the above parameters R, H{R) 
and /3o(-Rjj- According to Table 3, in the case of NP 0532 such values of the parameters /3o{R) and 
R correspond to the gravitational radiation power Lg w 4, 5 • 10^^ erg/sec. 
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Figure: 13. GMSW parameters in a neutron star magnetosphere 
for a dipole radial dependence of the magnetic field strength. The 
line corresponds to IgC^; the points correspond to IgT. 

As is seen from Fig. 13, the region favourable for the GMSW formation lies in the range 6i?4-16i?, 
i.e., where the local magnetic field strength is 3A10* 4A10^ G. When the GMSW pulse passes, 
these local values increase by a factor of 10 to 30. Thus a neutron star in whose magnetosphere 
a GMSW zone is formed, is able to radiate GW only from its magnetic poles, like pulsars with a 
pencil radiation pattern. In this case the probability of direct GW detection from such sources is 
drastically decreased. However, GMSW open another way for GW observation. A formed GMSW 
carries above all strong magnetic fields. They move from the neutron star in the magnetic equator 
plane, and therefore should lead to an increased pulsar magnetic bremsstrahlung intensity at the 
moment when the GMSW front passes. Thus anomalous electromagnetic radiation flashes in the 
pulsar radiation should be observed at moments when quadrupole oscillations are excited. 

The total magnetic bremsstrahlung intensity of a relativistic electron is proportional to squared 
magnetic field strength [lO] : 

I ^2e^H^f/i3m^c^), (83) 

where 

p = mv\J\ — V / 

is the electron momentum. 

Therefore the curves H^{S) shown in Figs. 9-12, actually describe the time dependence of the 
magnetosphere magnetic bremsstrahlung intensity, i.e. the local electromagnetic response to the 
gravitational radiation of the neutron star. Such a response might be detected by an observer at rest 
placed in the magnetosphere and screened from the electromagnetic radiation coming from other 
regions. The situation is more difficult with a total response of the magnetosphere to the GW, 
detected by a distant observer. We will later return to this problem. 

*In all further figures the magnetosphere is considered in the magnetic equator plane. 
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As was mentioned above, the response of a homogeneoiis magnetoactive plasma even to strietly 
periodic gravitational radiation has the form of a single pulse. But even if it were not the case, the 
response of a neutron star magnetosphere to a GW would still have the same form. Indeed, a shock 
wave (GMSW), emerging after the excitation of quadrupole oscillations of a neutron star, should 
throw the equatorial sector of the magnetosphere away into the interstellar space. For the next pulse 
to be formed, the magnetosphere should restore. The necessary time for its restoration is of the 
order of At ^ l/vs where I is the characteristic size of the magnetosphere and Vs is the velocity of 
sound. For a typical neutron star magnetosphere At ~ 1 sec. And typical quadrupole oscillation 
damping times comprise tenths of a second, according to Table 3. 



5.2 Effect of magnetosphere inhomogeneity upon GMSW 

According to Table 3, neutron star eigen-oscillation periods vary depending on the stellar mass in 
the range of 0,3 to 1,2 ms. Therefore the local duration of GMSW pulses should, by (79), satisfy the 
condition 

At < 7 • 10"^ ^ 3 • 10-*s, (84) (84) 

i.e., be shorter than 70 to 300 microseconds. Fig. 14 shows the dependence of the GMSW pulse local 
duration on the radial coordinate r, AT(r), calculated according to Eq. (78). 
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Figure: 14. The dependence of local pulse duration At (in mi- 
croseconds) on the distance form the star centre r/c (in microsec- 
onds). The duration was calculated by Eq. (78). The black 
circles mark the boundaries of a region where the GMSW effect 
is sufficiently well-developed. Outside this region the result is of 
a formal nature. On the left boundary of the region the pulse 
duration should quickly grow up to the values (84). 

As is seen from Fig. 14, the actual local pulse duration in the region where the GMSW mechanism 
is fairly effective, ranges from 1 to 10 /iS. Thus in this region At ~ 10~^r the following condition is 
fulfilled with a large spare: 



Atc <C r, 



(85) 
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justifying the use of the GMSW formulae for describing an inhomogeneous magnetosphere. 

An observer out of the neutron star magnetosphere would detect the magnetic bremsstrahlung 
from the magnetospheric electrons during the whole time while the local pulse passes through the 
magnctosphcric region r_ < r < r+ where favourable conditions for GMSW development arc realized, 
namely, (76), (77). With a certain caution these conditions may be specified: < 0, 5 (the lower 
bound of the GMSW range, r_ and T > 5 (its upper bound, r+). So the size of the GMSW zone is 

Ar = r+ — r_ . 

If Ar < 0, a GMSW zone does not appear in the neutron star magnetosphere at all. Since, as we 
have seen, the GMSW pulse spreading velocity is very close to c, the whole magnetic bremsstrahlung 
detected by remote observer will be concentrated in the time "window" of duration AT 

AT = Ar/c = t+-i_. (86) 

where t± = r± /c arc the instants when the GW leading front reaches the \ipper and lower boundaries 
of the GMSW zone. Near its boundaries r_ and r+ the GMSW is poorly developed (in the first 
case the first GMSW parameter is too large, in the second case the second parameter is too small). 
Therefore the intensity of the clcc;troniagnctic signal is small near the boundaries of the window, 
while in its medium domain a radiation maximum (large T and small ^^) is achieved. The form of 
the signal itself is yet to be calculated. Fig. 15 shows the dependence of the window width on the 
magnetic field strength and the GW magnitude. 
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Figure: 15. Dependence of the GMSW existence range AT (in 
fis ) in the magnetosphere of a neutron star of radius R= 1,2- 10® 
cm on the magnetic held strength H{R), (related to 10^^ G) and 
the GW magnitude (3o(R). The thin line corresponds to fio{R) = 
5 • 10"^ , the points to (io(R) = 7,5 • 10"^, the thick line to 
m = 10-^ 



5.3 Magnetic bremsstrahlung intensity 

Note that the local density of the bremsstrahlung intensity W(t,r) is determined by the local val- 
ues of the squared magnetic field strength, H^(t,r), and the local electron number density in the 
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magnetosphere, ne{t,r). For ultrarelativistic electrons by TIJ| it is 

where is the electron energy. Further we will assume that the size of a local pulse is much smaller 
than both the characteristic scale of the magnetospheric inhomogeneity r and the window width 
Ar. Thus the retarded time u is a quick variable and the radial coordinate r is a slow variable. 
Then a GMSW may be described by the formulae for homogeneous plasma, where it is necessary to 
use the local values of the GMSW parameters, ^^(r) and T(r). Mean while in the exact stationary 
solutions there arises a weak dependence on the radial coordinate r , i.e. the solution will be weakly 
nonstationary and the nonstationarity will show it selfin the form of a functional dependence of the 
GMSW solutions on the local values of the parameters, e.g., 

A(u;r) = 1 - T{r)q{u;r), 

etc. Thus, from the particle number conservation law and the solution stationarity it follows: 

L^n^ (r, t)v^ (u) = Const ~ -^n° (r) , (87) 

where n°{r) is the unperturbed electron number density in the magnetosphere. Taking into account 
that Eue = e/2 (half energy of a relativistic magnetosphere belongs to the electrons) and using the 
solutions (34), (35) and (45) for a weak GW, we obtain: 

n,(r,i) = n^(r)A-3/2(«); (88) 

E{r,t) ^ Eo{r)A-^^^u); (89) 
where EQ{r) is the unperturbed energy of magnetospheric electrons and 

H\u) ^ H^{r)A-^u). (90) 

Thus we get the following relation for the magnetic bremsstrahlung of the ultrarelativistic magneto- 
sphere: 

Wir,t) = Wo{r)A-^^/'^{u), (91) 

where Wo{r) is the magnetic bremsstrahlung intensity density for an nonperturbed magnetosphere. 
Eq. (91) needs a relativistic correction taking into account the plasma motion: the radiation density 
should be multiplied by the relativistic factor (1 — w^/c^)~^/^. The net result is 

W{r,t) = Wo(r)A^"/2l(A3/2 + A-3/2). (92) 

Thus W{r, t) ~ Wo{r)A^'^ . Note a large value of the exponent of the governing function A(w), which 
leads to a large steepness of the local magnetic bremsstrahlung pulse. Integrating Eq. (92) over 
the GMSW zone near the magnetic equator, we obtain a formula for the variation of the complete 
magnetospheric magnetic bremsstrahlung resulting from a GW pass: 

'■+ 

AJ{t) =271^0 J ^iA)Woirydr, (93) 
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where 

$(A) = iA-iV2(A3/2 + - 1; (94) 

©0 is the angle of the knife radiation pattern. This formula completely describes the shape of the 
signal to be detected by a remote observer. The expression in the square brackets in (93) is notably 
nonzero only in the domain of the local GMSW pulse, i.e. in the domain 

0<S ^TT-^ct- 7r/(wT(r)) <ct. (95) 

Therefore the integral (93) tends to zero for t < r^/c and t > r^/c + tt / (uiT (r)) . The observed 
pulse duration is formally determined by these limits. However, since, as noted above, near the zone 
boundaries the GMSW is poorly developed, the actually observed pulse duration (more precisely, its 
half-width) can turn out to be smaller than this value. Without solving the problem of the observed 
pulse shape, let us estimate its magnitude in its medium domain 

r-/c<t<r+/c, (96) 

when the pulse local duration is much smaller than the window width: 

At < AT. (97) 

Under these conditions the integrand in Eq. (93) is (5-like, therefore with a good precision the 
following estimate is valid: 

S- 

2ti < Hu(;-);-- > cB,, 



AJ{t) « "'^^ ^ J ^A{S))dS, (98) 

s+ 

where S± = ujT{t — r±/c), < Wo{r)r^ > and the value of r* is determined from the equation 
r* = ct- l/((wT(r*)). (99) 



Thus 



AJ(i) ~ 4^VFo(r*)r^cAT(A-^), (100) 
where 

(A-7) = A-'{S)dS; (101) 

s+ 

AS = ATloT. Thus, by order of magnitude, the intensity change of the magnetic bremsstrahlung 
(in its maximum) as a result of the GMSW excitation is equal to a product of the unperturbed 
magnetospheric radiation intensity in the GMSW zone by the dimensionless factor (A"'^) (which can 
reach 10~^). We will return to a calculation of the observed pulse shape in our next paper. Here we 
restrict ourselves to estimates of the complete energy of a GMSW pulse (Fig. 16). 
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Figure: 16. The radial coordinate (r) dependence of the pulse 
energy Epuise = ■inr'^Ss (in erg) hy Eq. (80). The GMSW effect is 
sufEciently well-developed only in the range 0, 8 < lg{r/R) < 1, 2 
(7, 6 • \{fcm < r < 1, 9 • lO^cmj. Outside this region the results 
are of clearly formal nature. The pulse energy should rapidly fall 
near the boundaries of this range. 



5.4 Source ofquadrup ole oscillations 

There naturally arises the question of a source ofpulsars' quadrupole oscillations. A possible energy 
source for such oscillations might be represented by explosive nuclear reactions with heavy hyperons 
like n + n < — > p + T,+ taking place in neutron stars cores at densities over 10^^ g/cm^ [13]. The pres- 
ence of strong magnetic fields should lead to an asymmetry of the explosions, i.e. to the quadrupole 
moment excitation. For such processes to occur in a neutron star, it must be sufficiently young. 
Numerical simulations of the process of a neutron star cooling shows [14] that after a Supernova 
explosion the neutron star temperature falls approximately by an order of magnitude in 10** years. 
Consequently, GMSW should be sought in radiation from sufficiently young pulsars formed no earlier 
than 10000 years ago. 



6 The Crab pulsar NP 0532 emits gravitational waves 

A pulsar with the required parameters does exist: it is the famous pulsar in the Crab nebula, NP 
0532, life-time is less than 1000 years (the 1054 Supernova). This pulsar is the youngest of all 
known ones (and consequently the hottest), it has the shortest period (at least among the closest 
pulsars, enumerated in Table 2): T = 0, 033 sec. What is suprising is that the radio emission of this 
pulsar contains anomalies which can be with a large degree of confidence identified with the GMSW. 
Namely: there are single irregular, the socalled giant pulses (on the average a pulse in every 5 to 10 
minutes) [11]. The radiation intensity in the giant pulses is a few tens of times higher (by roughly a 
factor of 60) than in common pulses. But the most interesting is that the duration of the giant pulses 
is no more than 9 • 10~^ s, i.e., almost by 2(!) orders shorter than that of the common pulses from 
NP 0532 (r ^ 6 • 10~^ s). The common pulse duration, as is easily seen, is about 1/7 of the NP 0532 
rotation period, so that the common pulses are clearly explained geometrically by pulsar rotation. 
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The giant pulse duration is 300 times shorter than the pulsar rotation period, and consequently the 
existence of the giant pulses has not yet found any satisfactory theoretical model. 

However, the giant pulses are easily explained by the GMSW, and their duration is not related 
to the pulsar rotation period, or an angle of the knife radiation pattern, but to its eigen-oscillation 
period, Tq. A comparison of the NP 0532 giant pulse duration with that of a GMSW pulse (84) 
shows a striking coincidence. Indeed, for the NP 0532 pulsar, as known from the annihilation line 
shift in the 7 radiation spectrum (400 keV instead of 511 keV), the gravitational redshift is known 

AE/E = MG/Rc^ = 0, 217. 

Then from Table 3 we find the pulsar mass: M = 1,67Mq and the corresponding neutron star 
radius: i? = 12 km. According to Eq. (79) and Table 3, the GMSW pulse duration for the NP 
0532 pulsar should be about 87 microseconds, while the observed NP 0532 giant pulse duration is 
approximately 90 microseconds ?? (!). 

For the Crab pulsar to — 2,5 ■ 10'^ years ??; then, setting SM — 1, 67, i? = 12 km, we find from 
(82): H ^ 5,2 ■ 10^^ G. The angle Qq of the knife radiation pattern is connected with the observed 
pulse duration r by the relation 89 = 2'kt/P. For the pulsar NP 0532 this angle is 0,48 « 28°. 
Assuming the complete pulsar luminosity in the continuous spectrum to be about 5 • 10^^ erg/s, we 
find the giant pulse intensity recalculated for the whole neutron star surface: 

Lg^ant ~ 4 • lO^^erg/s. 

As has been pointed out above, the real magnetic field strength on the pulsar NP 0532 surface is lO^^G 
|11] . To explain the observed giant pulse emission power, one needs GW magnitude values on the 
stellar surface of the order of 10~*. Note that, according to Fig. 15, the window width AT 180/^s 
corresponds observed GMSW pulse should be about 90 /is, which again precisely coincides with the 
observed giant pulse duration! 

The indicated magnitude I3q (R) corresponds to the gravitational radiation power of the order of 
4 • lO''^ erg/ s and the neutron star oscillations energy about E^ ~ 4 • 10^^ erg. In this case the neutron 
star surface oscillation magnitude is about 1 cm. Taking into account that in the whole lifetime of 
NP 0532 (1000 years) approximately 7- 10^ giant pulses have been emitted, we get an estimate of the 
energy carried away from the neutron star by GW for the whole time of its existence: i? = 2, 8 • 10^^ 
erg. It is 10~^ of the rest energy of this neutron star, which completely agrees with the assumption 
of a permanent rebuilding of its core. Thus, to a high degree of confidence we can state that the 
giant pulses observed in the pulsar NP 0532 radiation are optical manifestations of gravimagnetic 
shock waves (GMSW) excited by the gravitational radiation of the neutron star corresponding to 
the pulsar NP 0532 j20j. 

7 Conclusion 

Besides NP 0532, among all known pulsars only PSR 0833 seems to be able to emit (but more seldom) 
giant pulses. Other pulsars are too old for it. Therefore it is necessary to concentrate the main effort 
on observations of these two pulsars. It should be stressed that there is no other mechanism able 
to accelerate a shock wave to subluminal velocities. Therefore an investigation of the giant pulse 
spectrum in the X-ray range, aimed at discovering a violet shift in the radiation spectrum, is of 
utmost importance. A comprehensive study of the giant pulses (their shapes and instantaneous 
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spectrum) will allow one not only to verify the existence of gravitational radiation, but also to get 
additional information on the neutron stars structure and the processes in their interior. In turn it 
is necessary to study the GMSW pulse formation in detail theoretically. 
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